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LONG-TERM  GOALS 

A  key  question  in  space  weather  and  space  environment  study  is  to  understand  and  to  quantify  the 
contribution  to  thermospheric/ionospheric  variability  by  processes  from  both  the  lower  atmosphere  and 
from  the  magnetosphere  and  solar  input.  The  primary  goal  of  this  project  is  address  this  question  by 
developing  a  model  that  encompasses  the  entire  atmosphere,  from  the  ground/ocean  to  the  upper 
thermosphere.  We  will  achieve  this  goal  by  extending  the  NCAR  Whole  Atmosphere  Community 
Climate  Model  into  the  thermosphere  and  ionosphere  (WACCM-X),  making  use  of  the  physics  and 
many  of  the  algorithms  of  the  National  Center  for  Atmospheric  Research  Thermosphere-Ionosphere- 
Mesosphere-Electrodynamics  General-Circulation  Model  (NCAR  TIME-GCM).  Such  a  seamless 
whole  atmosphere  model  will  enable  us  to  self-consistently  study  how  elements  in  the  coupled  upper 
atmosphere/ionosphere  system  interact  with  one  another  and  to  determine  how  this  coupled  system 
responds  to  the  variable  energy  input  from  the  sun  and  the  complex  interactions  between  the  lower 
atmosphere/ocean  and  the  middle  and  upper  atmosphere.  The  information  from  this  research  will  be 
useful  for  ONR  to  develop  a  seamless  operational  model  that  simulates  the  present  day  structure  and 
dynamics  of  the  thermosphere-ionosphere-mesosphere-lower  atmosphere-system  including  its 
response  to  solar  variability  and  global  change. 

OBJECTIVES 

Our  scientific  objectives  are  to  study  the  nature  of  the  sources  of  variability  in  the  upper 
atmosphere/ionosphere  system  self-consistently  using  a  seamless  model  spanning  the  whole 
atmosphere.  The  primary  technical  objectives  of  the  project  are:  (1)  Extend  the  upward  boundary  of 
the  WACCM  to  the  top  of  the  thermosphere  (3x10'^  hPa,  around  ~500  km)  by  including  the  main 
thermospheric  and  ionospheric  processes.  (2)  Test  and  refine  the  WACCM-X  by  comparing  with 
empirical  models  and  with  observations  of  upper  atmospheric  tides,  ionospheric  electric  fields,  and 
geomagnetic  perturbations.  (3)  Help  the  development  of  NOGAPS-ALPHA  into  the  lower 
thermosphere  by  transferring  appropriate  packages  from  WACCM-X.  (4)  Help  the  development  of 
ionospheric  data  assimilation  schemes  at  the  Utah  State  University  by  providing  WACCM-X  output 
with  tidal  and  planetary  wave  variability. 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 
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APPROACH 


The  WACCM-X  is  built  upon  the  WACCM  version  3  (WACCM3),  which  is  based  on  the  NCAR 
Community  Atmospheric  Model  (CAM3)  and  now  a  component  of  the  Community  Climate  System 
Model  (CCSM),  and  includes  all  the  physical  parameterizations  of  that  model.  There  are  three 
dynamical  cores  available  in  CAM3,  and  the  finite  volume  dynamical  core  [Lin,  2004]  has  been  used 
for  WACCM3.  This  method  is  superior  for  mass  conservation  and  important  for  the  interactive 
chemistry  in  WACCM3.  The  WACCM3/WACCM-X  chemistry  module  is  derived  from  the  three- 
dimensional  chemical  transport  Model  for  Ozone  and  Related  chemical  Tracers  (MOZART)  [Brasseur 
et  ah,  1998],  which  solves  51  neutral  species  plus  5  major  ion  species  (02-I-,  0+,  NO+,  N+,  and  N2-I-) 
and  electron  density.  Processes  important  for  the  energetics  in  the  upper  atmosphere,  including 
photolysis  in  the  solar  spectral  range  from  the  EUV  to  750  nm,  photoionization,  solar  heating  at  short 
and  long  wavelengths,  including  infrared  transfer  under  non-local  thermodynamic  equilibrium  (NLTE) 
conditions,  auroral  processes,  ion  drag,  and  Joule  heating  are  calculated  or  parameterized  in 
WACCM3.  Gravity  wave  forcing  and  eddy  diffusion  induced  by  gravity  wave  breaking  are 
parameterized  in  WACCM3.  Molecular  diffusion  of  minor  species  diffusion  in  the  vertical  direction  is 
also  calculated.  A  detailed  summary  of  WACCM3  can  be  found  in  Garcia  et  al.  [2007]. 

The  present  upper  boundary  of  WACCM3,  at  around  ~150  km,  only  permits  studies  of  the 
thermosphere  up  to  about  -110-120  km,  because  higher  altitudes  are  significantly  affected  by  the 
assumed  upper  boundary  conditions.  The  first  phase  of  this  development  is  to  raise  the  upper  boundary 
of  the  model  to  the  upper  thermosphere  and  at  the  same  time  include  the  physics  that  is  relevant  for  the 
thermosphere.  In  particular,  major  species  diffusion  and  species  dependent  molecular  weight,  specific 
heats,  and  gas  constants  need  to  be  formulated  in  WACCM-X,  and  upper  boundary  conditions  for 
winds,  temperature,  and  various  species  appropriate  for  the  upper  thermosphere  need  to  be 
implemented.  With  these  developments,  WACCM  now  has  a  thermospheric  component.  The 
compositional,  thermal  and  wind  structures  are  being  validated  against  empirical  model  results  and 
TIME-GCM  results.  The  second  phase  of  the  development,  which  is  on-going,  is  to  develop 
ionospheric  electrodynamics  modules,  such  as  ambipolar  diffusion  and  ExB  drift,  in  WACCM-X. 

Also  funded  by  this  award,  we  have  further  developed  the  NCAR  thermosphere-ionosphere- 
mesosphere-electrodynamics  general  circulation  model  (TIME-GCM)  and  continued  to  use  the  model 
for  studying  the  upper  atmospheric  variability.  To  better  resolve  the  impact  of  lower  atmospheric 
variability  on  the  upper  atmosphere  and  to  take  advantage  of  the  augmented  computing  power  afforded 
by  the  new  NCAR  facilities,  the  spatial  resolution  of  TIME-GCM  has  been  doubled  in  each  dimension, 
and  the  lower  boundary  condition  of  TIME-GCM  can  now  be  specified  by  assimilated  data  from  the 
European  Centre  for  Medium-Range  Weather  Forecast  (ECMWF),  which  has  6-hourly  output.  This 
model  is  now  being  used  to  study  short-term  variability  due  to  planetary  waves  and  gravity  waves. 

NCAR  personnel  participating  in  this  work  include:  Han-Li  Liu  (dynamics  of  the  middle  and  upper 
atmosphere),  Raymond  Roble  (aeronomy  and  global  upper  atmosphere  and  ionospheric  dynamics), 
Arthur  Richmond  (electrodynamics  and  upper  atmosphere  waves),  Stanley  Solomon  (ionospheric 
physics),  Astrid  Maute  (scientific  support  on  ionospheric  physics),  Liying  Qian  (scientific  support  on 
thermospheric  physics),  Joseph  Mclnemey  (scientific/programming  support:  model  upward  extension), 
and  Benjamin  Foster  (programming  support  and  model  development).  WACCM  is  a  collaborative 
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effort  between  three  NCAR  Divisions:  The  High  Altitude  Observatory  (HAO),  The  Climate  and  Global 
Dynamics  Division  (CGD),  and  the  Atmospheric  Chemistry  Division  (ACD),  and  involves  scientists 
from  these  divisions  as  well  as  outside  collaborators. 

WORK  COMPLETED 

Extended  the  WACCM  (1.9°x2.5°  degrees  horizontal  resolution)  to  the  upper  thermosphere  at  3.4x10'^ 
hPa  (the  same  as  the  upper  boundary  of  the  TIME-GCM,  approximately  at  500  km).  In  the  process  of 
achieving  this  upward  extension,  we  have  (1)  Implemented  modules  to  resolve  the  major  species 
diffusion,  which  becomes  increasingly  important  above  110  km.  (2)  Implemented  modules  and  revised 
the  codes  to  reflect  the  constituent-dependency  of  the  specific  heats,  gas  constant,  and  mean  molecular 
weight.  They  were  set  to  constants  in  the  previous  codes.  (3)  Revised  the  treatment  of  the  vertical 
diffusion  equations  for  minor  species  and  heat  conduction  equation. 

Made  test  runs  using  the  extended  WACCM  over  one  model  year  under  solar  maximum,  solar  medium, 
and  solar  minimum  conditions,  and  obtained  the  wind  structure,  thermal  structure,  the  vertical  profiles 
of  major  species  (N2,  02,  O),  oxygen  species  (03,  OID,  02_1S,  02_1D),  hydrogen  species  (H,  H20, 
H02,  OH),  nitrogen  species  (NO,  N02),  carbon  species  (CO,  C02,  CH4),  ions  (0+,  02-I-,  NO+,  N+, 
N2+)  and  electron  density  from  the  ground  to  the  upper  thermosphere. 

Made  comparisons  with  TIME-GCM  and  MSIS/HWM  and  the  wind/temperature/compositional 
structures  in  the  middle  and  upper  atmosphere  are  in  good  agreement  with  those  obtained  from  TIME- 
GCM  and  MSIS/HWM.  This  is  the  first  time  that  the  wind,  temperature  and  compositional  structures 
are  resolved  from  the  ground  to  the  upper  thermosphere  self-consistently  in  a  seamless  model. 

Compared  the  thermospheric  density  under  solar  maximum,  medium,  and  minimum  conditions  from 
WACCM-X  simulations  with  MSIS.  WACCM-X  successfully  reproduced  the  semi-annual  variation  of 
the  0/N2  ratio,  with  both  the  magnitude  and  phase  in  good  agreement  with  MSIS. 

Compared  both  migrating  and  non-migrating  tides  from  WACCM-X  with  NASA  TIMED/SABER  and 
TIDI  observations.  The  seasonal  variation  of  these  tides  from  the  model  is  in  good  agreement  with  that 
from  the  observations.  The  thermospheric  tides  from  the  model  also  show  strong  short-term  variability. 

Migrated  WACCM3/WACCM-X  into  the  CCSM  framework,  so  that  the  models  can  now  run  with 
interactive  ocean  models  in  addition  to  specified  Sea  Surface  Temperature. 

Developed  and  tested  a  high  resolution  version  of  TIME-GCM  (2.5°x2.5°x0.25  scale  height)  with  the 
lower  boundary  specified  by  the  ECMWF  operational  model  results  at  10  hPa  (2.5°x2.5°  and  6-hourly 
output).  Initial  simulation  shows  that  the  tidal  variability  and  non-migration  tides  are  better  resolved. 

The  upper  boundary  of  the  TIE-GCM,  TIME-GCM  and  WACCM-X  are  at  approximately  500  km 
altitude.  To  self-consistently  model  the  ion  flux  and  ionospheric  electrodynamics  we  need  to  include  a 
model  of  the  plasmasphere.  For  this  purpose,  we  have  coupled  the  Global  Ionosphere  Plasmasphere 
Model  (GIP)  with  TIE-GCM  by  collaborating  with  George  Millward  from  the  NOAA  Space  Weather 
Prediction  Center.  The  modular  approach  in  coupling  the  TIE-GCM  and  GIP  makes  it  straightforward 
to  couple  WACCM  and  GIP  in  the  next  phase  of  development.  The  coupled  code  is  still  undergoing 
testing  of  all  the  interfaces  and  coupling  modes.  Immediate  future  work  will  include  comparing  the 
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results  of  the  coupled  GIP-TIEGCM  code  with  the  stand-alone  TIE-GCM  and  GIP  models  to  verify  the 
improvements. 

Further  development  of  TGCM/WACCM  postprocessors:  (1)  Further  developed  the  TGCM  IDE 
postprocessor  so  that  it  can  process  results  from  the  extended  WACCM  simulations.  (2)  Further 
developed  the  TGCM  FORTRAN  postprocessor  so  that  it  can  process  results  from  the  high  resolution 
version  of  TIME-GCM. 

Provided  the  Naval  Research  Laboratory  Navy  Operational  Global  Atmospheric  Prediction  System 
Advanced  Level  Physics-High  Altitude  (NRL  NOGAPS-ALPHA)  team  with  an  upgraded  version  of 
the  Global  Mean  Model  and  the  WACCM  FUV  module. 

Provided  the  Utah  State  University  Global  Assimilation  of  Ionospheric  Measurements  (USU  GAIM) 
team  with  WACCM  output  for  studying  ionospheric  data  assimilation  in  the  presence  of  variability  due 
to  lower  atmosphere  perturbations  (tides,  planetary  waves,  and  gravity  waves). 

Organized  a  CEDAR  workshop,  Ground-Space  Models  for  Studying  Atmospheric  Coupling  (Co- 
Chairs:  Han-Li  Liu,  Dave  Siskind  and  Bob  Schunk,  18  June,  2008,  Zermatt  Resort,  Midway,  Utah 
http://cedarweb.hao.ucar.edu/wiki/index.php/2008_Workshop:Ground- 

Space_Models_for_Studying_Atmospheric_Coupling).  The  two-session  workshop  was  intended  to 
review  the  current  status  of  such  models  and  to  discuss  model  development,  validation,  model  inter¬ 
comparison,  research  application  of  these  models,  and  integration  with  the  CEDAR  observational 
community. 

RESULTS 

Figure  1  shows  the  meridional  winds  at  March  equinox  from  WACCM-X  (top),  HWM  (middle),  and 
TIE-GCM  (bottom)  at  2.7x10'^  hPa  (~I20  km,  left)  and  6.8x10'^  hPa  (~330  km,  right)  under  solar 
maximum  conditions.  The  longitude  structures  and  the  diurnal  variation  of  the  meridional  wind  in  the 
upper  thermosphere  from  all  three  models  are  in  general  agreement,  and  the  wind  from  WACCM-X 
shows  more  variability.  At  lower  thermosphere,  winds  from  all  three  models  show  semi-diurnal  tidal 
structures  with  comparable  phase  and  amplitude.  In  addition,  the  WACCM-X  wind  has  a  clear 
wavenumber  3  (DE3)  signature  at  lower  latitudes,  indicating  the  presence  of  diurnal  eastward  3  non¬ 
migrating  tide  in  WACCM-X.  Figure  2  compares  the  DE3  temperature  and  wind  components  to  those 
measured  by  the  TIMED/SABER  [Garcia  and  Liebermann,  2007]  and  TIDI  [Wu  et  al.,  2008].  Both  the 
amplitude  and  seasonal  variation  of  these  components  compare  very  well  between  the  model  results 
and  the  observations. 

Figure  3  is  the  zonal  wind  in  the  altitude  range  between  the  ground  and  the  upper  thermosphere  from 
WACCM-X.  The  plot  clearly  shows  different  features  from  lower  to  upper  atmosphere:  below  the 
mesosphere  the  time  scale  is  dominated  by  planetary  waves,  in  the  upper  thermosphere  by  the 
migrating  tides,  and  in  between  the  temporal  and  spatial  scales  are  very  complex,  involving  planetary 
waves,  tides,  gravity  waves  and  interactions  of  these  components. 

Satellite  drag  data  show  that  there  are  annual  and  semi-annual  changes  of  thermospheric  density,  which 
is  not  well  resolved  in  WACCM  or  TIME-GCM.  Using  TIE-GCM,  Qian  et  al.  found  that  a  seasonal 
variation  of  eddy  diffusion  compatible  is  needed  to  obtain  the  density  variation  based  on  model 
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sensitivity  studies  and  satellite  drag  data.  Simulations  show  that  when  the  eddy  diffusion  with  seasonal 
variation  is  imposed  at  the  lower  boundary  of  the  TIE-GCM,  neutral  density  variation  consistent  with 
satellite  drag  data  and  column  0/N2  variation  consistent  with  measurements  by  TIMED/GUVI  are 
obtained.  These  model-data  comparisons  and  analyses  indicate  that  turbulent  mixing,  which  is 
primarily  related  to  gravity  waves  from  the  lower  atmosphere,  may  contribute  to  seasonal  variation  in 
the  thermosphere,  particularly  the  asymmetry  between  solstices  that  cannot  be  explained  by  other 
mechanisms.  (Qian,  L.,  S.  C.  Solomon,  and  T.  J.  Kane,  Seasonal  variation  of  thermospheric  density  in 
response  to  forcing  by  the  lower  atmosphere,  J.  Geophys.  Res.,  submitted,  2008.) 

Liu  et  al.  used  WACOM  to  examine  the  chaotic  error  growth  in  the  whole  atmosphere  context,  which 
will  help  us  better  understand  the  predictability  of  the  whole  atmosphere  system  and  future  application 
of  data  assimilation  to  forecasting  the  system.  From  ensemble  WACCM  simulations,  they  found  that 
the  early  growth  of  differences  in  initial  conditions  are  associated  with  gravity  waves,  this  becomes 
apparent  first  in  the  upper  atmosphere  and  progresses  downward.  The  differences  later  become  more 
profound  on  increasingly  larger  scales,  and  the  growth  rates  of  the  differences  change  in  various 
atmospheric  regions  and  with  seasons,  and  correspond  closely  with  the  strength  of  planetary  waves. 
The  growth  rates,  on  the  other  hand,  are  not  sensitive  to  the  altitude  where  the  small  differences  are 
introduced  in  the  initial  conditions  or  the  physical  nature  of  the  differences.  Furthermore,  the  growth 
rates  in  the  middle  and  upper  atmosphere  are  significantly  reduced  if  the  lower  atmosphere  is  regularly 
reinitialized,  and  the  reduction  depends  on  the  frequency  and  the  altitude  range  of  the  re-initialization. 
(Liu,  H.-L.,  F.  Sassi,  and  R.  R.  Garcia,  Error  growth  in  a  whole  atmosphere  climate  model,  J.  Atmos. 
Sci.,  in  press,  2008.) 

Deng  et  al.  examined  the  non-hydrostatic  effect  on  the  upper  atmosphere  by  comparing  results  from 
the  Global  Ionosphere  and  Thermosphere  Model  (GITM)  and  TIME-GCM.  Their  results  show  that 
after  a  sudden  intense  enhancement  of  high-latitude  Joule  heating,  the  vertical  pressure  gradient  force 
can  locally  be  25%  larger  than  gravity,  resulting  in  a  significant  disturbance  away  from  hydrostatic 
equilibrium.  This  disturbance  is  transported  from  the  lower  altitude  source  region  to  high  altitudes  by 
an  acoustic  wave,  which  has  been  simulated  in  a  global  circulation  model  for  the  first  time.  Due  to  the 
conservation  of  perturbation  energy,  the  magnitude  of  the  vertical  wind  perturbation  increases  with 
altitude  and  reaches  150  (250)  m/s  at  300  (430)  km  during  the  disturbance.  The  upward  neutral  wind 
lifts  the  atmosphere  and  raises  the  neutral  density  at  high  altitudes  by  more  than  100%.  These  large 
vertical  winds  are  not  typically  reproduced  by  hydrostatic  models  of  the  thermosphere  and  ionosphere. 
The  results  give  an  explanation  to  the  cause  of  such  strong  vertical  winds  reported  in  previous 
observations.  (Deng,  Y.,  A.  D.  Richmond,  A.  J.  Ridley,  and  H.-L.  Liu,  Assessment  of  the  non¬ 
hydrostatic  effect  on  the  upper  atmosphere  using  a  general  circulation  model  (GCM),  Geophys.  Res. 
Lett.,  35,  L01104,  10.1029/2007GL032182,  2008.) 

Using  simulation  results  from  WACCM,  Liu  et  al.  revisited  the  wind  balance  in  the  mesosphere  and 
lower  thermosphere  and  demonstrated  that  the  geostrophic  balance  approximate  holds  in  the 
meridional  direction  but  is  no  longer  valid  in  the  zonal  direction  due  to  the  large  zonal  gravity  wave 
forcing.  As  a  result,  the  zonal  mean  geostrophic  meridional  wind  is  significantly  different  from  the 
actual  zonal  mean  meridional  wind,  and  the  residual  mean  merdional  circulation  derived  from 
geostrophic  winds  is  much  weaker  than  that  derived  from  model  winds.  The  ageostrophic  contribution 
comes  primarily  from  gravity  wave  forcing,  so  that  there  is  an  approximate  three-way  balance  between 
pressure  gradient,  Coriolis  force,  and  gravity  wave  forcing.  The  relationship  between  the  geostrophic 
winds  and  actual  winds  is  also  tested  using  measurements  from  TIMED  SABER  and  TIDI  instruments. 
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This  research  provides  a  possible  route  to  directly  infer  gravity  wave  forcing  from  meridional  wind 
measurements.  (H.-L.  Liu,  D.  R.  Marsh,  Q.  Wu,  and  J.  Xu,  Wind  balance  in  the  mesosphere  and  lower 
thermosphere,  in  preparation.) 

Yuan  et  al.  used  multiple  years  of  continuous  observations  of  mesopause  temperature  and  horizontal 
wind,  each  lasting  longer  than  24  hours,  from  Colorado  State  University  Na  Lidar  Facility  Fort  Collins, 
CO  (41°N,  105°W)  to  derive  monthly  climatology  and  seasonal  variations  of  temperature  and 
horizontal  winds.  The  observed  mean-state  in  temperature,  zonal  and  meridional  winds  are  compared 
with  the  predictions  of  3  current  general  circulation  models:  WACCM3,  the  Hamburg  Model  of  the 
Neutral  and  Ionized  Atmosphere  (HAMMONIA),  and  the  year  2003  simulation  of  the  TIME-GCM. 
While  general  agreement  is  found  between  observation  and  model  predictions,  there  exist 
discrepancies  between  model  prediction  and  observation,  as  well  as  among  predictions  from  different 
models.  The  altitude  of  winter  zonal  wind  reversal  and  seasonal  asymmetry  of  the  pole-to-pole 
meridional  flow  are  also  compared,  with  the  importance  of  full-diumal-cycle  observation  for  the 
determination  of  mean  states  studied.  (T.  Yuan,  C.-Y.  She,  D.  A.  Krueger,  F.  Sassi,  R.  R.  Garcia,  R.  G. 
Roble,  H.-L.  Liu,  and  H.  Schmidt,  Climatology  of  mesopause  region  temperature,  zonal  wind  and 
meridional  wind  over  Fort  Collins,  CO  (41°N,  105“W)  and  comparison  with  model  simulations,  J. 
Geophys.  Res.,\\2>,  D03105,  2008.) 

Marsh  et  al.  used  WACCM3  to  study  the  atmospheric  response  from  the  surface  to  the  lower 
thermosphere  to  changes  in  solar  and  geomagnetic  forcing  over  the  1 1 -year  solar  cycle.  Energy  inputs 
include  solar  radiation  and  energetic  particles,  which  vary  significantly  over  the  solar  cycle.  This  paper 
presents  a  comparison  of  simulations  for  solar  cycle  maximum  and  solar  cycle  minimum  conditions. 
Changes  in  composition  and  dynamical  variables  are  clearly  seen  in  the  middle  and  upper  atmosphere, 
and  these  in  turn  affect  the  terms  in  the  energy  budget.  Generally  good  agreement  is  found  between  the 
model  response  and  that  derived  from  satellite  observations.  A  small  but  statistically  significant 
response  is  predicted  in  tropospheric  winds  and  temperatures  which  is  consistent  with  signal  observed 
in  reanalysis  data  sets.  (Marsh,  D.  R.,  R.  R.  Garcia,  D.  E.  Kinnison,  B.  A.  Boville,  K.  Matthes,  F.  Sassi, 
and  S.  C.  Solomon.  Modeling  the  whole  atmosphere  response  to  solar  cycle  changes  in  radiative  and 
geomagnetic  forcing.  J.  Geophys.  Res.,  112,  D23306,  2007.) 

Liu  studied  the  large  wind  shear  and  fast  transport  above  the  mesopause  and  possible  relationship  with 
gravity  waves  and  tides.  With  the  maximum  atmospheric  static  stability  above  the  mesopause,  large 
wind  shears  can  be  obtained  under  the  limit  of  dynamic  stability  and  the  consistently  large  observed 
wind  shears  there  can  thus  be  related  to  the  stability  constraint  set  by  the  background  atmosphere.  The 
maximum  wind  shears  determined  from  this  stability  limit  are  found  to  be  in  general  agreement  with 
the  outer  envelope  of  most  of  the  large  wind  shears  from  chemical  release  experiments  at  low  and  mid¬ 
latitudes.  Diagnostic  calculations  also  indicate  that  the  meridional  transport  in  this  region  may  not  be 
well  understood  solely  by  examining  the  mean  meridional  circulation,  and  large  amplitude 
tides/planetary  waves  can  play  an  important  role  in  the  bulk  transport  of  tracers.  Strong  stochastic 
winds,  presumably  due  to  gravity  waves,  do  not  seem  to  significantly  change  the  large  scale  pattern  of 
the  transport  but  may  extend  the  range  of  the  tracer  movement.  (Liu,  H.-L.,  On  the  large  wind  shear 
and  fast  meridional  transport  above  the  mesopause,  Geophys.  Res.  Lett.,  34,  L08815, 
doi:10.1029/2006GL028789,  2007.) 

Liu  et  al.  studied  the  large  mesospheric  temperature  inversion  and  coincident  strong  short-term  tidal 
variability  using  temperature  and  wind  measurements  from  lidar,  NASA  TIMED/SABER  and  TIDI 
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instruments  and  compare  these  results  with  TIME-GCM.  With  a  large  transient  planetary  wave 
specified  at  the  model  lower  boundary,  the  model  is  able  to  produce  strong  diurnal  tidal  variability 
comparable  to  that  from  the  lidar  observation,  and  the  modeled  temperature  inversion  is  similar  to  that 
from  the  SABER  measurement.  The  model  results  suggest  that  the  planetary/tidal  wave  interaction 
excites  non-migrating  tides  and  modulates  the  gravity  modes  and/or  the  rotational  modes  of  the  diurnal 
migrating  tide.  Among  the  non-migrating  tides,  the  diurnal  zonally  symmetric  (S=0)  component  is  the 
strongest,  and  its  interaction  with  the  planetary  wave  leads  to  a  strong  diurnal  eastward  wavenumber  1 
component.  (Liu,  H.-L.,  T.  Li,  C.-Y.  She,  J.  Oberheide,  Q.  Wu,  M.  E.  Hagan,  J.  Xu,  R.  G.  Roble,  M.  G. 
Mlynczak,  J.  M.  Russell  III,  Comparative  study  of  short  term  tidal  variability,  J.  Geophys.  Res.,  112, 
doi:  10.1 029/2007 JD008542,  2007.) 

IMPACT/APPLICATIONS 

WACCM,  TIE-GCM,  and  TIME-GCM  are  community  models  and  the  development,  evaluation  and 
application  of  these  models  have  and  still  will  have  extensive  community  involvement.  The  WACCM3 
has  been  recently  released  (http://waccm.acd.ucar.edu/)  and  can  be  downloaded  from  the  NCAR  data 
portal  (http://cdp.ucar.edu/).  The  TIE-GCM  will  soon  be  released  through  the  Community  Coordinated 
Modeling  Center  (CCMC).  In  addition,  we  participate  in  NRL  studies,  NSF  Coupling  and  Energetics 
of  Atmospheric  Regions  (CEDAR),  and  Space  Weather  Initiative  (SWI)  programs.  We  worked  closely 
with  NASA  TIMED  team  to  help  interpret  the  observational  results,  and  a  NASA  Small  Explorer  team 
for  mission  planning.  The  WACCM  provides  a  seamless  modeling  tool  to  study  whole  atmosphere 
coupling  and  connecting  the  tropospheric  climate  with  space  weather,  and  it  attracts  participation  from 
the  lower,  middle  and  upper  atmosphere  communities.  We  also  participate  in  the  NCAR  Climate 
Systems  Modeling  effort  in  examining  the  couplings  between  the  upper  and  lower  atmospheres  and  in 
an  attempt  to  understand  the  effects  of  the  variable  solar  outputs  on  the  coupled  Earth  system.  We  also 
work  closely  with  other  seamless  modeling  teams,  such  at  the  NRL  NOGAPS-ALPHA  team. 

RELATED  PROJECTS 

The  numerical  modeling  effort  is  complemented  by  a  data  analysis  and  interpretation  efforts,  and 
comparative  studies  with  observations  and  other  model  development.  We  have  actively  participated  or 
have  participation  from  the  following  projects/missions: 

NASA  TIMED  Satellite  Mission. 

NASA  Living  With  a  Star  Program. 

NASA  Geospace  SR&T  Program. 

NSF  CEDAR:  Whole- Atmosphere  Modeling  of  the  Thermosphere/Ionosphere  Responses  to  Lower- 
Atmosphere  Dynamics  and  Variability. 

NSF  Lidar  Consortium. 

NCAR  Community  Climate  System  Model. 

NRL  NOGAPS-ALPHA. 

USU  GAIM. 
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Figure  1:  Meridional  wind  at  March  equinox  front  WACCM-X  (top),  HWM  (middle) 
and  TIE-GCM  (bottom)  at  2.  hPa  (~120  km,  left) 

and  6.8x10'  hPa  (~330  km,  right)  under  solar  maximum  conditions. 
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Figure  2:  Temperature  (top)  and  zonal  wind  (bottom)  components  of  diurnal  eastward  wavenumber 
3  (DE3)  from  WACCM-X (left)  and  TIMED/SABER  (top  right)  and  TIDI  (bottom  right). 

The  pressure  heights  of  the  model  results  have  been  chosen  to  be  close  to  those  of  SABER  and  TIDI. 
The  TIDI  plot  has  been  rotated  to  make  the  layout  the  same  as  the  other  plots. 
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Figure  3:  Zonal  wind  from  WACCM-X from  the  ground  to 
the  upper  thermosphere  at  52N  on  December  23. 
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